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conduct  a hardware  requirements  analysis  and  trade-off.  In  particular,  the 
paper  addresses  the  application  of  the  general  framework  for  an  aircraft  to 
satellite  precision  tracking  problem. 
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Introduction 


There  are  maiy  aerospace  problems  which  include  the  requirement  for 
precision  pointing  and  tracking  from  one  accelerating  vehicle  to  another. 

If  knowledge  of  the  position  and  velocity  states  of  the  two  vehicles  is 
perfect,  the  pointing  and  tracking  problem  is  straightforward.  Unfortunately, 
such  ljiowledge  is  never  available.  In  practice,  t'ue  instrumentation  on 
boar<\  the  tracking  vehicle  is  imperfect  and  provides  measurement  information 
which  is  corrupted  by  noise  and  other  unwanted  effects.  For  example,  a 
radar  tracking  device  will  in  general  be  inperfect  because  of  random  phase 
front  distortion  or  scintillation  noise  which  might  be  target  dependent  and 
could  prevent  the  required  degree  of  accuracy  from  being  achieved.  In 
addition,  the  target  vehicle  dynamics  might  be  unknown  and  must  therefore 
be  modeled  as  a random  process.  For  example,  when  passively  tracking 
another  aircraft,  the  aircraft  acceleration  profile  is  unknown  but  could  be 
modeled  as  a first  order  Markov  process  [9,10,11  J.  The  model  for  a low 
orbit  satellite  could  include  the  low  order  gravitational  harmonics  for  the 
liarth;  high  order  harmonics  and  other  perturbations  could  then  be  modeled 
as  random  processes.  In  general,  therefore,  perfect  knowledge  of  the  system 
state  is  not  available  in  the  pointing  and  tracking  pioblem,  and  stochastic 
estimation  techniques  are  necessary. 

This  paper  contains  the  development  of  a general  framework  for  applica- 
tion of  the  extended  Kalman  filter  to  the  precision  pointing  and  tracking 
problem  in  order  to  estimate  certain  necessary  physical  variables  in  the 
problem.  The  fundamental  system  dynamics  are  described  by  a set  of  primary 
state  equations  to  which  are  adjoined  the  state  equations  necessary  to 
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describe  instrument  errors.  This  total  state  description  is  the  "true 
system?*  model.  However  this  model  is  of  high  order  (typically  60-100 
states)  and  the  optima*  filter  based  upon  it  could  not  therefore  be  iaple- 
mented  with  a small  on-board  computer.  On  the  other  hand,  the  true  system 
model  or  truth  model  may  be  used  in  a covariance  sensitivity  analysis  to 
yield  the  true  performance  of  a reduced  order  filter  which  is  small  enough 
for  on-board  implementation.  By  testing  the  sensitivity  of  approximations 
from  the  truth  model  to  the  reduced  order  filter,  a hardware  trade-off 
analysis  will  indicate  instrumentation  requirements.  The  reason  for  this 
appro  at.  i along  with  the  necessary  sensitivity  equations  is  more  fully 
developed  in  reference  [21].  This  general  approach  is  applied  in  the  paper 
to  the*  particular  problem  in  which  a hi$i  altitude  aircraft  must  track  a 
low  polar  orbit  satellite  to  precise  degree  of  accuracy.  Some  of  the 
results  showing  the  sensitivity  of  the  tracking  errors  to  measuring  instru- 
ment precision  are  presented.  More  generally,  the  framework  has  been 
applied  by  the  authors  or  could  be  applied  r.  aircraft  to  aircraft,  moving 
grouid  vehicle  tc  aircraft,  and  satellite  to  satellite  problems. 

Work  in  the  area  of  Kalman  filtering  for  pointing  and  tracking  includes 
that  by  Fitts  [1,2, 3, 4]  in  which  filtering  is  accomplished  in  the  inertial 
reference  frame.  Detailed  modeling  of  error  sources  was  not  considered. 

A general  study  by  Pearson  [5,6]  again  does  not  include  error  source  modeling. 
Other  studies  in  this  area  by  Landau  [7],  Fitzgerald  [8,9]  and  others 
[10,11,12,13]  do  not  include  error  source  modeling.  None  of  the  references 
considers  the  particular  problem  of  aircraft  to  satellite  tracking. 
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General  Framework 
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The  fundamental  objective  is  to  establish  an  accurate  line  of  sight  to 
the  target,  vehicle  being  tracked,  along  which  some  fora  of  electromagnetic 
link  will  be  maintained.  For  example,  to  illuainate  a target  with  a very 
narrow  heamwidth  ccwraxii cations  laser  could  require  tracking  accuracies  in 
the  order  of  microradians.  Even  with  highly  sophisticated  measuring  and 
control  instruments  it  may  not  be  possible  to  achieve  sudi  accuracy  without 
using  optimal  estimation  and  control. 

This  section  considers  the  modeling  of  the  physical  variables  necessary 
to  establish  a line  of  sight  in  both  position  and  velocity  and  the  modeling 
of  the  instrument  error  sources  inherent  in  the  problem  of  estimation  for 
pointing  and  tracking.  The  approach  is  general  as  it  may  be  used  for  a 
number  of  applications  ot  which  only  one  is  shown  in  the  next  section. 

In  this  general  approach  it  ir  assumed  that  the  target  is  passive. 

That  is,  it  does  not  assist  the  tracker  in  any  way  by  providing  target  state 
information.  This  assumption  implies  that  some  model  must  be  developed  for 
the  target  to  obtain  the  estimate  of  target  acceleration  which  is  necessary 
to  solve  the  tracking  problem. 

Figure  1 illustrates  the  basic  tracker  geometry.  The  tracking  device 
is  controlled  in  azimuth  and  elevation  relative  to  a reference  coordinate 
system.  The  tracker  coordinate  system  is  shown  in  this  figure.  Note  that 
the  xt  axis  aligns  with  the  tracking  device  bore-sight.  For  perfect 
tracking  this  will  lie  along  the  line  of  sight  between  tracker  and  target. 
Figure  2 illustrates  the  two  coordinate  frames  of  primary  interest.  The 
tracker  frame  is  misaligned  from  the  line  of  sight  frame  because  of  inperfect 
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tracking.  Note  that  this  Misalignment  is  described  by  two  Euler  «igle  rota- 
tions 5c  and  6n  which  are  chosen  in  general  to  be  about  axes  where  Mis- 
alignment measurements  would  in  practice  be  available  from  a sensor. 

The  need  for  a target  acceleration  estimate  to  solve  the  tracking 
problem  has  already  been  indicated.  To  control  the  tracking  device  it  is 
necessary  to  estimate  the  line  of  sight  angular  velocity  and  the  tracker 
misalignment.  These  estimates  provide  rate  and  position  feedback  control 
respectively  to  the  tracker  motors.  In  some  problems  it  may  also  be  neces- 
sary to  estimate  range  and  range  rate  of  the  target.  Now  with  the  exception 
of  target  acceleration,  all  these  parameters  are  line  of  sight  and  tracker 
parameters  and  since  ideally  the  tracker  frame  will  coincide  with  the  line 
of  sight  frame  it  is  logical  to  dioose  the  latter  fr»s*  in  which  to  model 
the  system.  This  choice  has  the  disadvantage  that  target  acceleration  must 
be  traisformed  from  inertial  coordinates  in  which  estimates  are  most  likely 
to  be  available,  into  tracker  or  line  of  sight  coordinates.  An  alternative 
would  be  to  model  completely  in  the  inertial  frame  but  this  would  necessitate 
transforming  the  line  of  sight  and  tracker  states  into  inertial  coordinates. 
This  paper  uses  the  line  of  sight  frame  for  modeling  since  this  involves 
the  least  nimfccr  of  coordinate  transformations. 

The  system  dynamics  is  singly  described  by  considering  the  relative 
position  vector  of  the  target  from  the  tracker.  Let  r be  the  relative 
position  vector  and  the  angular  velocity  vector  of  the  line  of  sight 

relative  to  inertial  space.  Then,  by  the  Theorem  of  Cariolis, 
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where  indicates  the  differential  is  taken  with  respect  to  inertial 
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Finally,  the  motion  between  the  tracker  coordinate  frame  and  the  line  of 
sight  coordinate  frame  is  characterized  by  the  following  matrix  differential 
equation  (151: 
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where  C 5 is  the  coordinate  transformation  matrix  front  lino  of  sight 
coordinates  to  tracker  coordinates  and  Wjg,  Wt  are  the  line  of  sight  and 
tracker  an^.lar  velocity  skew  symmetric  "cross  product”  matrices  respectively, 
liquations  (1),  (2)  and  (3)  are  used  to  obtain  the  primary  system  state  equa- 
tions excluding  target  state  equations  as  follows: 

From  equation  (1)  let  R = range  * Irj  and  let  V = range  rate,  then: 

R - Vr  (4) 


Let  ar  he  the  relative  acceleration  vector  of  target  from  tracker  which  is 
equal  to  * ^RACKFJT  “*  hc  tJie  anPt,lar  velocity  relative  to 

inertial  space,  subscripts  x,  v,  and  z indicate  vector  components  in 
a right  hand  orthogonal  coordinate  system,  superscripts  Is  and  t indicate 
the  relevant  coordinate  system,  and  the  small  angle  approximations  to  be 
used  as 
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From  equations  (2)  and  (3): 
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Equations  (4-9)  are  the  state  equations  describing  the  range,  range 

rate,  lino  of  sight  angular  velocity  and  tracking  device  misalignment.  Note 

Is 

that  the  component  of  lines  of  sight  angular  velocity  is  not  included 

sx 

as  a state,  since  wgular  velocity  along  the  line  of  sight  is  not  required 
for  control.  However  the  tracker  angular  velocity  vector  (u^,  which  is 
included  as  a system  parameter,  is  required  along  all  three  tracker  axes. 

The  relative  acceleration  vector  a^  is  the  difference  between  the  target 
acceleration  vector  atar  and  the  tracker  acceleration  vector  a^  which  is 
a system  parameter.  The  state  modeling  for  target  acceleration  will  be 
described  below.  It  is  inportant  to  note  at  this  point  that  the  inclusion 
of  system  parameters  such  as  the  angular  rate  of  the  tracker  implies  that 
measurements  of  those  parameters  are  available. 

The  modeling  of  the  target  acceleration  is  the  most  difficult  part  of 
the  problem.  The  reader  will  recall  that  no  target  state  information  is 
directly  available  since  the  target  is  assumed  generally  to  be  passive.  If 
the  target  vehicle  is  an  aircraft,  then  an  exponentially  time  correlated 
acceleration  model,  which  implies  that  the  target  motion  is  uncertain,  may 
be  used  [9,llj.  Clearly,  I r more  information  about  the  target  acceleration 
is  available,  then  this  should  he  used  in  the  model.  I'or  example,  if  it  is 
certain  that  the  target  does  not  maneuver,  then  a zero  acceleration  model 
is  appropriate. 

If  the  target  is  a satellite,  then  accurate  equations  of  motion  may  be 
developed  and  utilized  to  give  an  acceleration  model.  A.,  initial  orbit  may 
be  given  through  ephemcris  data.  This  might  be  no  more  than  an  initial  orbit 
prediction  necessary  for  satellite  acquisition.  The  gravitational  motion 
due  to  the  non- spherical  Earth,  the  Siai  and  the  Moon  can  be  accurately 
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modeled  although  the  model  may  be  complex.  The  effects  of  drag  and  solar 
radiation  pressure  are  more  difficult  to  describe  L.nce  they  depend  on  the 
physical  chaiacteri sties  of  the  satellite,  i.e.,  its  mass,  surface  area 
and  shape.  Furthermore,  if  the  satellite  is  unknown  then  there  may  be  a 
finite  probability  that  the  stellite  is  independently  comnanded  to  execute 
a AV  maneuver  or  attitude  change.  It  may,  in  fact,  be  acceptable  to 
ignore  all  these  effects  but  they  must  initially  be  included  in  the  truth 
model  until  the  filter  sensitivity  to  them  has  been  deteimined  for  a 
particular  orbital  condition. 

Perfect  knowledge  of  the  satellite  physical  characteristics  cannot  be 
assumed  in  the  general  approach  to  the  problem.  It  is  reasonable  however 
to  assume  some  statistical  knowledge.  To  describe  the  atmospheric  drag  and 
solar  radiation  pressure  effects,  the  use  of  a ballistic  coefficient  and 
solar  pressure  coeffic?nt  is  necessary.  Since  the  satellite  physical 
characteristics  are  inlikel>  to  change  with  time  these  coefficients  can 
reasonably  be  modeled  as  random  biases  with  the  general  state  equation 

X = 0. 

Satellite  AV  maneuvers  will  most  probably  be  of  short  duration 
unless,  for  exanple,  the  satellite  is  continually  thrustin^  to  follow  a 
drag- free  orbit.  Modeling  is  best  accomplished  based  on  some  a priori 
knowledge  of  the  satellite  function  and  adaptive  techniques.  Alternatively, 
a simpler  approach  is  to  assun?  a worst  case  size  and  probability  of 
maneuver  and  thereby  retain  a large  measure  of  uncertainty  in  the  orbit. 

This  safe  approach  can  however  unnecessarily  degrade  the  ultimate  tracking 
capability  if  the  maneuver  does  not  occur. 
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In  sunnary,  each  type  of  target  vehicle  presents  a different  problem 
the  physics  of  which  must  be  carefully  analyzed.  If  no  information  is 
available  then  the  exponentially  correlated  randcm  variable  model  [9,11] 
may  oe  used.  Any  information  which  is  available  must  be  used  if  the  best 
model  is  to  be  found. 

The  problem  of  measurement  errors  was  briefly  described  in  the  intro- 
duction but  the  extent  of  this  problem  was  not  explained.  When  a reduced 
order  filter  is  applied  in  the  real  world  it:;  performance  can  be  highly 
sensitive  to  measurement  errors.  To  inplement  a reduced  order  filter  it  is 
therefore  necessary  to  conduct  a full  sensitivity  analysis  which  implies 
that  the  system  truth  model  must  include  all  known  error  sources.  The  sbb- 
scquerc  analysis  will  hopefully  show  that  maiy  of  the  errors  can  be  ignored 
or  discarded  and  what  ine  expected  filter  degradation  will  be. 

Particular  applications  of  this  general  framework  generate  slightly 
different  measurement  requirements.  The  following  if.  a description  of  a 
typical  measurement  set  which  will  suit  most  high  precision  requirements. 
Measurements  are  assumed  to  be  available  for  the  following  states  and 
parameters: 

tracker  angular  velocity,  w*; 

tracker  origin  inertial  acceleration,  at; 

tracking  device  misalignment  from  the  1 rue  line  of  sight,  6e,  Sn; 
range,  R;  and 
range  rate , V . 

The  length  of  this  paper  prohibits  a detailed  description  of  every  measure- 
ment, but  tn  illustrate  the  principles,  consider  the  measurement  model  of 
tracker  velocity.  The  three  rate  gyroscopes  are  mounted  along  the  three 


tracker  axes  and  yield  the  Measurement 


o».  * w + 6k  u + b + c 

X*  L • L • . • 

1 m i gi  1 gi  gi 


i * x,y,z. 

Now  the  components  of  are  cotx,  and  wtz  which  are  the  true 
angular  velocity  components  along  the  ti acker  x,  y and  z axes.  The 
measurement  vector  components  are  u>txjn,  and  u>tzm  respectively. 

The  terms  <5k^  are  torque  scale  factor  errors.  During  manufacture  and 
assembly  of  the  rate  gyro,  the  calibration  process  will  remove  the  torque 
scale  factor  error  as  far  as  possible  but  there  will  remain  a small  residual 
component  which  cannot  be  compensated  for.  Since  such  a scale  factor  error 
is  unlikely  to  change  with  time,  particularly  in  the  short  term,  a random 
bias  model  is  chosen.  This  model  has  the  differential  equation  for  the 
x-gyro  scale  factor  error  for  example  of 


0. 


The  statistics  of  this  factor  may  be  found  from  gyro  test  data  or  from  an 
engineering  estimate.  Note  that  the  random  bias  mode  inplies  that  both 
mean  l:|(SkgX|  and  second  moment  E J <5^-gX2 { stay  constant  with  time. 

The  drift  characteristic  of  a rate  gyro  can  be  expressed  as  the  sum  of 
a constant  bias  drift  and  a time  varying  drift.  The  terms  b . , i * x,  y,  z 
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are  bias  drifts  which  describe  taose  components  of  gyro  drift  which  do  not 
change  in  time.  The  time  varying  drifts  are  represented  by  the  terms 
c^,  i * x,  y,  z.  Again,  the  bias  terms  b ^ are  modeled  using  the  bias 
random  variable,  for  exanple 


V * 


The  model  used  for  the  c^-  terms  would  depend  on  the  type  and  quality 
of  gyro.  Typically  such  time  correlated  drifts  can  be  described  by  expo- 
nentially time  correlated  random  variables  with  the  first  order  differential 
equation  form,  for  example 


c * - S c + v2  8 o u 

gx  pgx  gx  v gx  gx  gx 


in  which  Bnv  * and  t is  the  drift  process  correlation  time,  a 

gx  t^x  gx  gx 

is  the  rms  value  of  the  process  and  u isaunity  variance  white  driving 

noise.  Note  again  however  that  the  model  requires  knowledge  < f the  drift 

process  in  terms  of  process  correlation  time  and  rms  value. 

A rate  gyro  is  ideally  insensitive  to  accelerations.  However,  in  prac- 
tice, a mass  unbalance  will  exist  because  the  center  of  mass  and  center  of 
rotation  do  not  coincide.  The  mass  unbalance  coefficients  k-^  do,  however, 
remain  constant  with  time  and  can  therefore  again  be  represented  with  the 
random  bias  model , for  exarrple 


k = 0 
xy 


where  k is  the  coefficient  of  mass  unbalance  along  the  tracker 
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x- direction  due  to  the  mass  unbalance  effect  along  the  tracker  y- direction. 

If  the  rate  gyroscope  is  also  likely  to  be  sensitive  to  products  of  accelera- 
tion or  g2-sensitive  errors  then  the  model  should  also  include  random  bias 
coefficients  to  account  for  these  effects. 

In  practice,  the  axes  along  which  the  rate  gyroscopes  are  sensitive  do 
not  coincide  with  the  true  tracker  coordinates.  A measurement  from  one  rate 
gyro  will  therefore  include  conponents,  however  small,  of  angular  velocity 
along  the  other  two  tracker  axes.  By  assuming  that  the  displacement  r gles 

are  small  enough  for  small  angle  approximations  to  be  valid,  the  coordinate 

T 

trais  format  ion  matrix  can  be  approximated  as  follows.  is  the  coordinate 
transformation  matrix  from  gyro  coordinates  (G)  to  trackei  coordinat  ) 
where 


where 


Thus  tnc  term 


T 1 

A w.  , l = 

L J 1 
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accounts 

formation  and  since  misalignment  will  be  constant  with  time  the  terms 

can  be  described  by  the  random  bias  model,  for  example 


0. 


Finally,  an  additive  white  noise  £ j , i * x,  y,  z is  included  in  the  model 


to  account  for  those  higher  order  effects  which  are  not  otheiwise  modeled. 
Choice  of  variance  for  this  noise  may  be  difficult  to  determine  but  should 
be  of  the  same  order  of  magnitude  as  that  of  the  smallest  modeled  effect. 

At  this  point  it  is  worthwhile  reminding  the  reader  that  the  above 
revue  1 includes  all  reasonable  and  common  sources  of  error.  A sensitivity 
analysis  will  indicate  the  dominant  error  terms,  but  note  that  each  problem 
and  associated  measurement  sensor  will  produce  different  results.  Hopefully, 
the  effect  of  many  of  the  error  terms  in  a particular  application  will  be 
small  relative  to  the  required  accuracy  and  only  the  dominant  errors  will 
therefore  need  to  be  included  in  the  filter. 

The  following  is  a summary  of  the  remaining  measurement  models  with  a 
brief  description  of  each  term. 


Tracker  Acceleration.  The  measurement  model  for  tracker  acceleration  is 


m 


ati  * 6kai  ati  + bai  + cai  + kil  ati 

+ ki2  ati  + ki3  atj  + ki4  atk 

+ A di  a + s . ; i = x,y,z  ; j ,k^i 
L A 111  31  ;jj*k 


(ID 


in  which  at.  i = x,  v,  z are  the  measurements  of  the  true  tracker  origin 
lm 

acceleration  components  a^,  i * x,  y,  z.  The  accelerometer  scale  factor 
errors  k . and  bias  errors  b . are  modeled  again  as  random  biases.  c_; 
are  accelerometer  drift  errors  modeled  as  first  order  exponentially  corre- 
lated random  variables.  The  coefficients  k^j,  k^*  k.j  and  k^  are 
non-linear  calibration  coefficients  to  account  for  gJ,  g3  and  cross  track 


13 


accelerations  respectively  and  are  modeled  as  random  biases.  is  an 
error  angle  ttans format  ion  matrix  relating  acceleration  in  actual  (A) 
accelerometer  coordinates  to  acceleration  in  nominal  (N)  accelerometer 
coordinates.  The  elements  of  this  matrix  are  small  constant  angles  modeled 
as  random  biases.  In  general,  the  nominal  accelerometer  coordinate  system 
does  not  align  with  the  tracker  coordinate  system.  If  necessary,  the  trans- 
formation 


■ ^ 


T _N 


can  be  accomplished.  Finally,  is  an  additive  white  noise  term  to 

account  for  the  remaining  unmodeled  errors. 


Tracker  Angular  Misalignment  6c,  6n  (Bores ight  Error) 
The  tracker  misalignment  measurement  model  is 


k (6  + S)  + 5k  + b + £ 

e e f r.  e e 

k (6  ♦ s ) ♦ 6k  + b ♦ i 
nv  n n n n n 


(12) 

(13) 


where  k and  k are  deterministic  scale  factors,  S and  S are 

e n c n 

target  induced  scintillation  errors  as  observed  from  the  tracker  frame. 

These  errors  can  be  modeled  as  coupled  first  order  exponentially  time- 
correlated  random  variables  with  equations 

(14) 

(15) 


K = ' Sc  Se  * wtx  Sn  + °z  uc 


S = - 8 S - w S ♦ y/2  8 a u 

n nn  txevnnn 
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The  5k  and  6k  scale  factor  errors  and  b raid  t-  bias  errors  are 
e n c 

modeled  using  the  simple  bias  random  variable  model.  Finally,  t and  l 

t*  * i 

are  additive  white  noise  terms  to  account  for  unmodeled  effects. 

Range.  The  range  measurement  model  is 

Rm  = kr  (R  + Sr5  + kr  R + br  + Cr  (16) 

where  kr  is  a deterministic  scale  factor,  Sr  is  a range  scintillation 
error  modeled  as  a first  order  exponentially  correlated  random  variable, 
k , a scale  factor  error  and  br,  a bias  error  are  modeled  as  random 
biases.  Sj.  is  an  additive  white  noise  to  account  for  unmodeled  effects. 

Rajige  Ratc.  The  range  rate  measurement  model  is 

Vr  & kv  (Vr  + V + kv  Vr  + \ + «v4  (17> 

m 

where  the  terms  have  appropriate  significance  as  for  range. 

The  additional  state  equations  resulting  from  the  above  measurement 
models  arc  secondary  state  equations  which  are  adjoined  to  the  primary 
state  description  to  form  the  full  state  truth  model.  For  the  measurements 
described  there  are  63  secondary  state  equations.  For  the  typical  tracking 
problem  therefore  the  total  state  dimension  for  the  truth  model  is  about  75 
depending  upon  the  number  of  states  required  to  model  the  target  acceleration. 
While  it  is  doubtful  that  any  practical  problem  would  require  the  inclusion 
of  all  the  above  secondary  state  equations  in  the  truth  model,  even  a modest 
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requirement  to  adjoin  secondary  states  results  in  a large  dimension  truth 
model  which  could  not  be  handled  in  an  on-board  ccnputational  facility. 

It  is  inportant  to  examine  each  measurement  and  decide  hew  that  informa- 
tion should  be  used.  Clearly,  the  measurements  of  6£,  6^,  R and  Vr 

are  direct  measurements  of  primary  system  states  and  will  therefore  form  a 
system  measurement  vector.  The  measurement  01  a^  however  is  a parameter 
value  which  is  substituted  into  the  appropriate  state  equations.  The  measure- 
ment of  u)t  could  he  handled  in  an  identical  way.  However,  in  this  case  the 
information  can  be  considered  as  a pseudo-measurement  of  the  angular  velocity 
of  the  line  of  sight.  The  justification  for  this  reasoning  is  that  the  track- 
ing device  will  ideally  follow  the  line  of  sight  in  such  a way  that  the  mean 
angular  velocity  deviation  between  and  is  zero. 

So  far,  the  system  truth  model  state  and  measurement  equations  have  been 


developed.  Two  problems  are  apparent;  the  equations  arc  in  general  nonlinear 
ami  the  vector  dimensions  are  too  great  for  the  optimal  filter  to  be  imple- 
mented on  board  the  tracking  vehicle.  The  nonlinearity  of  the  equations  is 
overcome  by  linearizing  so  that  the  equations  for  the  extended  Kalman  filter 
are  valid.  It  is  necessary  however  to  implement  a filter  on  board  the  track- 
ing vehicle.  A reduced  order  (dimension)  filter  must  therefore  be  designed. 
In  gener.il  the  performance  of  a reduced  order  filter  will  be  sub-optimal 
and  must  be  evaluated  against  the  theoretically  optimal  performance  of  a 
filter  based  on  the  full  truth  state  description.  This  evaluation  is  carried 
out  initially  using  the  method  of  covariance  sensitivity  analysis  described 
by  Asher  and  Reeves  [19 | . However,  it  can  also  be  shown  [20]  that  the 
reduced  order  filter  is  in  general  conditionally  biased  and  to  obtain  all 
the  statistical  information  the  results  of  reference  [20]  nwst  be  applied. 
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The  process  of  designing  a reduced  order  filter  using  the  covariance 
sensitivity  analysis  will  not  only  indicate  which  errors  and  perturbations 
can  be  ignored  but  will  also  show  the  degree  of  precision  necessary  to 
achieve  a specified  tracking  accuracy  with  a particular  reduced  order 
filter,  there  is  currently  no  mathematically  precise  algorithm  by  which 
such  a design  process  can  be  conducted,  but  the  flow  chart  in  Figure  3 
shows  the  basic  method  applied.  The  filter  should  also  be  simulated  via 
Monte  Carlo  methods. 

An  initial  intuitive  choice  of  reduced  order  filter  parameters  and 
states  is  made.  For  example,  the  gyro  measurement  might  be  assumed  to  take 
the  simple  form 


tim 


= Cl) 


ti 


♦ c 


Ri 


+ r . 
’gl 


where  the  predominant  drift  term  has  been  retained  hut  the  remaining  terms 
have  been  accounted  for  by  increasing  the  variance  of  the  additional  white 
noise  £ . . For  another  example,  a satellite  orbit  may  be  assumed  to  be  a 
simple  two  body  orbit  in  which  the  gravitational  harmonics  and  other  pertur- 
bations are  modeled  by  a simple  additive  white  noise  of  suitable  variance. 
With  this  basis,  the  covariance  sensitivity  analysis  is  used  to  tune  the 
filter  by  adjusting  the  reduced  order  filter  parameters  so  that  the  estima- 
tion error  and  filter  sensitivity  are  minimized.  Note  that  this  will  be  an 
iterative  procedure  but  no  changes  are  made  to  the  system  truth  model  at 
this  time. 

The  tuning  process  will  also  indicate  where  the  filter  is  most  sensitive 
to  measurement  error  sources.  An  examination  of  this  sensitivity  is  then 
made  and  areas  for  possible  hardware  changes  or  trade-offs  are  identified. 
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For  example,  it  may  be  apparent  that  specifications  on  gyro  drift  could  be 
relaxed  but  angle  track  scinti'  lation  errors  must  be  reduced.  The  system 
is  therefore  redesigned  and  changes  are  made  to  both  the  full  state  optimal 
filter  aid  the  reduced  order  sub-optimal  filter  before  the  tuning  process 
is  repeated. 


Application  in  Aircraft -to- Satellite  Tracking 


To  illustrate  the  application  of  the  general  framework  described  above, 
consider  the  problem  in  which  a high  altitude  aircraft  is  tracking  a low 
polar-orbit  satellite.  Assume  the  required  accuracy  is  arbitrarily  taken 
to  be  on  the  order  of  microradiais.  Applying  a covariance  sensitivity 
analysis  will  show  first,  which  states  must  be  retained  in  a reduced  order 
filter,  and  secondly,  what  the  hardware  requirements  might  be  in  order  to 
meet  or  improve  on  the  accuracy  requirement. 

The  system  truth  model  comprises  the  state  and  measurement  equations 
described  in  the  general  framework  together  with  a set  of  target  state 
equations  which  describe  the  motion  of  the  satellite.  The  target  equations 


are 
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B * 0 (20) 

S - 0 (21) 


where  xg  and  vg  are  the  satellite  inertial  position  and  velocity  vectors. 
The  gravitational  force  due  to  the  Earth  is  described  by  the  vector  a^ 
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and  for  this  problem  includes  ha 'monies  up  to  sixth  order.  The  atmospheric 
drag  force  is  described  by  the  vector  a^  which  is  dependent  on  a satellite 
ballistic  coefficient  B.  The  perturbational  forces  due  to  the  sun  and 
moon  are  described  by  the  vectors  ag  and  am  respectively  and  the  solar 
pressure  force  is  described  by  the  vector  ap  which  is  dependent  on  a solar 
pressure  coefficient  S.  Finally,  the  uimodeled  effects  of  higher  order 
gravitational  harmonics  and  other  perturbations  are  accounted  for  by  the 
additive  noise  vector  which  is  assuiied  to  be  a zero  mean  white  noise. 
Assuming  that  the  satellite  is  passive  (but  not  cooperative),  the  probability 
that  the  satellite  executes  any  AV  maneuver  or  attitude  adjustment  is 
zero.  The  physical  properties  of  the  satellite  such  as  mass  shape,  size 
and  surface  area  are  assumed  to  remain  constant  for  the  duration  of  tracking. 
Thus  the  best  physical  knowledge  available  for  the  satellite  is  assumed  to 
be  a mean  and  standard  deviation  to  describe  an  expected  distribution  of 
these  satellite  physical  properties.  With  this  basis,  the  two  coefficients 
K and  S are  chosen  to  be  nuidom  biases  with  the  state  description  given 
by  equations  (20)  and  (21). 

The  resulting  truth  model  for  this  application  has  a state  vector  dimen- 
sion of  61,  a measurement  vector  dimension  of  5 and  a parameter  measurement 
vector  for  tracker  angular  velocity  and  tracker  origin  acceleration  of 
dimension  6.  The  accelerometer  crosstrack  errors  were  assumed  small  enough 
to  be  neglected.  The  accelerometer  and  gyro  bias  terms  were  included  into 
the  initial  conditions  for  the  respective  correlated  drift  terms.  A reduced 
order  filter  is  chosen  by  making  the  following  simplifications  from  the 
truth  model. 
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a.  Satellite  motion  can  be  described  by  a simple  two -body  orbit  per- 
turbed by  a white  d>ivM  ? noise  vector.  The  variance  of  this  vector  (£s 
in  the  true."  r xlel)  is  increased  to  accouit  for  the  two-body  approximations. 

h.  The  equations  describing  line  of  sight  angular  velocity  and  range 
rate  are  simplified  to  remove  the  terms  (see  equations  5,  6 and  9)  which 
account  for  the  effect  of  tracker  boresight  error.  In  each  equation,  a 
white  driving  noise  of  appropriate  variance  is  introduced  to  account  for 
the  approximation. 

c.  In  each  measurement  equation,  all  sources  of  error  are  removed  and 
the  variance  of  the  white  driving  noise  is  increased  in  each  case  to  compen- 
sate for  the  approximation. 

The  following  set  of  state  and  measurement  equations  results  from  the 
first  attempt  at  a reduced  order  filter. 


State  liquations 
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Measurement  Equations 
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The  objective  of  tuning  a reduced  order  filter  using  the  covariance 
sensitivity  analysis  technique  is  to  force  the  reduced  order  filter  error 
covariance  to  track  the  true  full  state  filter  error  covariance  as  closely 
as  possible.  In  practice  this  is  achieved  by  adjusting  the  various  filter 
noise  parameters  while  maintaining  the  truth  model  parameters  fixed  and  is 
a lang  process.  If  the  reduced  order  filter  has  been  over- sinplified  and 
significant  error  sources  are  not  modeled,  then  the  filter  may  be  extremely 
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difficult  to  tine  or  more  likely  will  show  a divergent  performance 
characteristic. 

Consider  then  the  filter  described  above.  For  the  particular  orbit 
chosep  and  the  initial  measuring  instrument  parameter  set,  several  changes 
are  necessary  before  the  filter  can  be  timed.  Because  of  the  low  orbit 
profile,  atmospheric  drag  produces  a significant  perturbation  from  the 
two-body  trajectory  and  cannot  be  properly  accouited  for  with  white  noise. 

The  drag  vector  a^  imrst  therefore  be  re-introduced  although  it  is 
sufficient  to  assume  a fixed  value  for  the  ballistic  coefficient  B.  It 
is  also  necessary  to  re-introduce  the  two  states  describing  target  - induced 
angle  track  scintillation  to  obtain  acceptable  reduced  order  filter  perfor- 
mance. In  this  case,  the  target- induced  angle  track  scintillation  noise  is 
assumed  to  have  a 2-second  correlation  time.  Because  this  error  has  a 
significant  steady  sthte  standard  deviation  compared  *o  other  error  sources 
in  the  angle  track  measurement,  it  is  not  possible  to  account  for  it  in  the 
additive  white  noise  ♦erm.  The  covariance  analysis  for  the  initial  choice 
of  reduced  order  filter  therefore  results  ir  an  increased  filter  state 
dimension  from  12  states  to  14. 

To  better  illustrate  this  tuning  process,  see  Figures  3 and  4.  Figure  3 
shows  the  standard  deviations  (1  - sigma)  of  the  tracker  misalignment  angle 
Se  erro1'-  predicted  by  the  reduced  order  filter.  Ine  filter  apparently 
performs  well  and  a steady  state  standard  deviation  of  approximately  10  yrad 
results.  Without  carrying  out  a covariance  analysis  the  engineer  might 
consider  this  to  be  very  satisfactory.  However,  Figure  4 is  the  result  of 
the  covariance  analysis  and  shows  the  true  standard  deviation  of 
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committed  by  the  particular  choice  of  reduced  order  filter  which  produced  j 

Figure  3.  In  fact,  as  already  explained,  it  is  necessary  to  re-introduce  j 

the  angle  track  scintillation  error  in  the  reduced  order  filter  before  I 

> 

satisfactory  performance  can  be  achieved.  Figures  5 and  6 on  the  other 

hand  illustrate  much  better  tuning  and  in  this  case  range  error  standard  j ' 

deviation  is  shown.  The  filter  is  slightly  over-estimating  the  true  error.  j 

I \ 

The  next  stage  in  the  sensitivity  analysis  is  to  find  the  error  sources  j 

included  in  the  truth  model  which  have  the  most  significant  effect  on  track- 
ing accuracy.  For  this  problem  and  an  assumed  parameter  set,  it  is  clear 
from  results  that  rate  gyro  drift  and  angle  track  scintillation  are  the  two 
predominant  error  sources.  To  illustrate  this  insight,  consider  the  remain-  ] 

ing  measurements  first.  A significant  error  standard  deviation  in  range 
measurement  scintillation,  bias  and  noise  can  he  tolerated  before  angle 
track  accuracy  begins’  to  suffer.  Range  rate  measurement  is  fouid  to  be  4 

super fluo> is  provided  range  information  is  available  at  2 sec  or  smaller 
intervals.  The  accelerometers  can  be  of  average  precision  capable  of  j 

measuring  to  within  about  1”  of  true  acceleration.  As  rate  gyro  drift  is 
improved  however,  the  tracking  error  standard  deviations  inprove  considerably,  i 

but  sensitivity  diminishes  until  improvements  beyond  an  error  standard  devia- 
tion of  about  0.5  * 10  B rad/sec  have  little  effec  . In  fact,  the  rate  gyros 

j 

could  theoretically  be  perfect  and  no  further  inprovement  would  result  \ 

j 

because  the  tracking  accuracy  is  ultimately  limited  by  angle  track  scintilla-  , 

tion  noise.  Figure  7 shows  a family  of  sensitivity  curves  for  tracking  1 

i 

accuracy  against  gyro  drift  and  angle  track  scintillation. 
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This  analysis  is  based  ipon  an  assuned  nominal  trajectory  which  is 
used  to  evaluate  the  system  matrices.  This  is  an  approximate  procedure 
that  is  valid  for  small  errors  and  is  exact  only  for  the  case  of  linear 
dynamics.  However,  the  procedure  is  extensively  and  successfully  used  for 
filter  analysis. 

The  problem  is  clearly  conplex  and  cannot  be  solved  by  intuition  alone. 
In  the  above  simplified  description  of  a sensitivity  analysis  such  improve- 
ment in  an  error  source  implies  a hardware  change.  The  hardware  change 
correspondingly  results  in  a truth  model  parameter  change  which  carries  the 
penalty  of  re-tuning  or  re-designing  the  reduced  order  filter.  Moreover,  as 
the  sensitivity  analysis  proceeds  and  truth  model  adjustments  (with  hardware 
implications)  are  made,  the  sensitivity  characteristic  changes.  For  example, 
if  in  this  problem  the  rate  gyro  and  angle  track  measurements  are  both 
perfect,  the  tracking  -accuracy  improves  but  becomes  sensitive  to  other 
measurements  aid  states  such  as  range,  acceleiation  and  the  satellite 
orbital  estimate.  An  error  budget  may  be  used  to  show  the  relative  benefit 
in  overall  performance  of  changing  one  sensor  versus  another,  so  that  cost 
effective  hardware  decisions  may  be  made. 

Figure  8 shows  a flow  chart  for  a typical  filter  performance  evaluation. 
The  inner  loop  illustrates  the  reduced  order  filter  tuning  and  re-designing 
to  achieve  satisfactory  performance  against  a specific  truth  model.  The 
outer  loop  illustrates  the  hardware  requirement  and  trade-off  process  through 
which  the  system  truth  model  is  adjusted  to  reflect  changes  in  hardware. 
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Conclusions 


A general  framework  has  been  developed  for  the  application  of  estimation 
techniques  to  the  precision  pointing  and  tracking  problem  from  one  accelera- 
ting vehicle  to  another.  Modeling  is  carried  out  in  the  line  of  sight 
coordinate  system  in  which  measurement  information  is  most  likely  to  be 
available.  There  is  no  general  model  for  target  motion  but  it  can  be  assumed 
that  target  motion  will  most  easily  be  described  in  inertial  coordinates. 
Specific  problems  with  differing  targets  will  require  different  target  models. 

The  problem  of  handling  large  computational  requirements  has  been 
identified  and  hence  the  need  to  find  a reduced  order  filter  and  tune  this 
filter  against  a full  state  system  truth  model.  Furthermore,  the  method  of 
identifying  hardware  requirements  and  trade-offs  to  meet  specific  performance 
criteria  using  the  covariance  sensitivity  analysis  has  been  described. 

Finally,  the  general  framework  has  been  applied  to  the  specific  problem 
of  tracking  a satellite  from  an  aircraft  to  an  arbitrary  degree  of  precision. 
The  particular  problem  involved  a high  altitude  aircraft  and  a low  polar 
orhit  satellite.  Some  results  following  from  the  covariance  sensitivity 
analysis  for  this  problem  have  been  presented. 
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Appendix  A 


Gravitational  Potential  Model 

The  model  for  the  gravitational  potential  is 


where  kg  is  the  gravitational  constant  for  the  earth,  m is  the  mass  of 
the  earth,  r is  the  radial  distance  of  the  body  from  the  earth  center, 

are  Legendre  functions,  A^  is  the  longitude  of  the  satellite  with 
respect  to  the  Greenwich  medium,  and  m and  m are  harmonic  coef- 
ficients for  the  potential  model. 


Drag  Force 

The  model  for  the  drag  is 

Ad 

where  p is  the  atmospheric  density,  assumed  exponential,  B is  the  vehicle 
ballistic  coefficients,  V„  is  the  magnitude  of  vehicle  velocity  relative 

a 

to  the  rotating  atmosphere,  and  r is  the  velocity  vector  of  the  vehicle 

o 

relative  to  the  rotating  atmosphere. 


In  B 


Solar  Pressure 

The  model  for  the  solar  pressure  is 
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- KS 


whore  (•)  is  the  coordinate  of  the  sun  relative  to  the  vehicle,  K 

VS 

a proportionality  constant,  is  the  distance  from  the  sun  to  the 

vehicle,  and  S is  the  solar  pressure  coefficient. 
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TIME 


Fig  3*  Reduced  order  filter  predicted  performance  - 6 
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Fie  4.  True  reduced  order  filter  performance  - <$e 
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Fig  5.  Reduced  order  filter  predicted  performance  - Rai 
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Fig  6.  True  reduced  order  filter  performance  - Range 


